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INTRODUCTION 

Previous Technical Reports covered the period up to 14 October, 1991. In this Final 
Report, we summarize the principle accomplishments of this project. Appendix 1 is 
a paper submitted for publication in the Special AGU Monograph on results from the 
NASA CDP Program. The research described in this report is supported in part by 
NSF and NASA/LAGEOS II. 


SUMMARY 

The long-term goal of this project is to undertake a detailed study of the 
consequences of the Arabian plate convergence against Eurasia and its effects on 
the tectonics of Anatolia and surrounding regions of the eastern Mediterranean. A 
primary source of information is time rates of change of baseline lengths and relative 
heights determined by repeated SLR measurements (Sellers and Rands, 1990; 
Nooman et al., 1990; Ellmer et al., 1990). These SLR observations are augmented 
by a network of GPS stations in Anatolia, Aegea, and Greece, established and 
twice surveyed since 1988. The existing SLR and GPS networks provide the spatial 
resolution necessary to reveal the details of ongoing tectonic processes in this area 
of continental collision. Our effort has involved examining the state of stress in the 
lithosphere and relative plate motions as revealed by these space-based geodetic 
measurements, seismicity, and earthquake mechanisms as well as the aseismic 
deformations of the plates from conventional geodetic data and geological evidence. 
These observations are used to constrain theoretical calculations of the relative 
effects of: (1) the push of the Arabian plate; (2) high topography of Eastern Anatolia 
(-2000 m); (3) the geometry and properties of the African-Eurasian plate boundary; 
(4) subduction under the Hellenic Arc and southwestern Turkey; and (5) internal 
deformation and rotation of the Anatolian plate. 

ACCOMPLISHMENTS 

Our principal efforts have been directed along the following lines: (1) Continuing 
analysis of seismicity and earthquake focal mechanisms; (2) Expanding our finite 
element models to include the entire Anatolian plate, the Aegean Sea and the 
Northeastern Mediterranean Sea, northwest Iran and the Caucasus; 

(3) Investigating the tectonic processes associated with ongoing continental 
collision in eastern Anatolia using neo-tectonic and geophysical information; and 

(4) Initiating an international program to densify SLR coverage in the eastern 
Mediterranean with GPS (our participation in this effort has been partly supported 
by NSF and NASA Lageos II). These efforts, described below, are designed to 
provide an improved basis for interpreting NASA’s SLR measurements in this 
region. 



Seismicity. Earthquake Source Mechanisms, and Conventional Geodetic 
Observations 


The seismicity of the Middle East and that of Turkey are shown in Figures 2 and 4, 
respectively. Epicenters generally outline the major fault zones shown in Figures 1 
and 3. However, there are a significant number of earthquakes occurring away from 
these “plate boundaries”, indicating the presence of smaller faults and deformations 
in the plate interiors. 

Earthquake source mechanisms (Toksoz et al., 1977; McKenzie, 1978; 
Buyukasikoglu, 1979, 1980; Kocaefe and Ataman, 1982; Rotstein and Kafka, 1982; 
Jackson and McKenzie, 1984; Lyberis, 1984; Eyidogan and Jackson, 1985; Alptekin 
et al., 1986; England and Jackson, 1990) shown in Figure 5 define the relative 
motions across the faults and provide information on prevailing stress directions. 

The fault mechanisms are consistent with “tensions” in the Aegean region. This 
indicates a major change in the state of stress in the Anatolian plate. We are 
continuing to compile and re-analyze earthquake mechanisms as well as the 
geologic evidence in central Anatolia to determine where and how the stress regime 
changes. Mechanisms for extension in Aegea are being coupled with other 
mechanisms to explain this change in the state of stress. One possibility to be 
taken into account is the rotation of the Anatolian plate. Paleomagnetic data from 
Central Anatolia (Sanver and Ponat, 1981) show that the Anatolian plate has 
rotated counter clockwise by about 15° in the past 45 million years. Such a rotation, 
coupled with a possible anchoring due to subduction could contribute to the tensile 
regime of Western Anatolia. Westaway (1990) shows evidence for the rotation of 
individual blocks in western Turkey. These models are being incorporated as 
constraints in both close scrutiny of the data and theoretical modeling of plate 
interactions (see below). 

In addition to deformations shown by earthquakes, there is ample evidence for 
aseismic deformation in and around the Anatolian plate. Geodetic and creep 
measurements on the North Anatolian fault show that there is creep on at least one 
segment of the fault, at a rate of about 1 cm/yr. To the west of the site, a major 
geodetic network shows that one branch of the fault (that broke in the 1967 
earthquake) is locked, while a southern branch is deforming at a rate of 1.2 cm/year. 
The combination of seismic dislocations and creep gives a fault dislocation rate of 
about 1.2 to 2.0 cm/year across the North Anatolian Fault zone. This deformation 
rate is in good agreement with rates derived from our analysis of SLR and GPS 
observations (see below and Appendix 1). These are the first hard numbers we 
have on the present-day motion of the Anatolian plate and are extremely important 
for constraining models of contemporary plate interactions. 

Finite Element Modeling 

The purpose of this aspect of our research program is to develop a simple but robust 
approach to modeling plate interactions using a contact problem in 2-D elastic finite 
element method. An individual plate is considered as a continuum, whereas an 



aggregate of plates is treated as a discontinuum such that plate boundaries are 
represented as contact surfaces. The behavior of plates at the contacts is defined 
by the Coulomb-Navier failure criterion, and three types of contacts are considered: 
sliding contact (slip), tension release (separation), and sticking contact (single- 
node continuum). The first and second modes correspond to transcurrent and 
divergent plate boundaries respectively. Convergent motion at plate boundaries is 
achieved by the double-node differential displacement technique with a recursive 
solution. Preliminary models for the Eastern Mediterranean, where tectonic 
deformation is produced by the interactions of the Eurasian, Arabian and African 
plates, was carried out. Figure 6 shows the finite element grid and defined contact 
zones used for these preliminary models. GPS estimates of horizontal crustal 
motions across the North Anatolian fault and in western Turkey (see below and 
Appendix 1) and geological observations of fault offsets were used to constrain 
displacements. Figure 7 shows modeled displacements for the Eastern 
Mediterranean region. We find that the deformation pattern in the Eastern 
Mediterranean is substantially controlled by differential motion between the African 
and Arabian plates and southwestward migration of the Hellenic arc. We note that 
the southwestward migration of the Hellenic arc required to fit the GPS 
observations is roughly consistent in sense with initial estimates of crustal motion 
determined from SLR observations (Ambrosius, personal communication 1991). 

The larger motions of the arc suggested by these models (relative to those 
estimated from SLR) may reflect SLR and/or GPS measurement uncertainties and/or 
the influence of buoyancy forces on crustal motions in western Turkey. Additional 
modeling is being undertaken to determine the relative importance of buoyancy 
forces caused by elevation changes. 

Neotectonic Studies 


1. Strike-Slip Fault Geometry and Its Relationship to Earthquake Activity (Barka et 
al., 1987, 1988; Barka and Kadinsky-Cade, 1988). 

Analysis of Turkish strike-slip faults indicates that detailed fault geometry plays an 
important role in controlling earthquake rupture. Empirical relationships are used to 
estimate possible locations and sizes of future strike-slip events. These results 
have implications for earthquake activity on other strike-slip faults such as the San 
Andreas in California. 

2. Tectonic Processes of the Continental Collision in the Vicinity 
of the Maras Triple Junction, Southern Turkey (Barka et al., 1990) 

Study of the geometry, kinematics and slip rates of the structures surrounding the 
Maras Triple Junction, formed by the collision between the Arabian, African and 
Eurasian plates in eastern Turkey, provides information on the continental collision 
process. Deformation is diffused over a region more than 500 km wide, from the 
Taurus to Caucasus Mountains within which several continental blocks escape from 
the maximum convergence zone. The collision initiated about 13 Ma ago and had 
two phases: the first caused thickening of the crust by thrusting, the second 



involved block escape. We estimate that 40 % of the overall shortening is taken up 
by block escape along strike-slip faults, and the remainder by thrusting and folding 
throughout the collision zone. 

3. Eastern Turkey and Caucasia (Alptekin et al., 1986; Kadinsky-Cade et al., 1989) 

Armenian (7 December 1988) and Iran (20 June 1990) earthquake source 
mechanisms provide good examples of the deformation patterns of the lesser 
Caucasus and the western shores of the Caspian Sea, respectively. The 
mechanisms (large “beach balls” in Figure 8) indicate that both earthquakes 
included a combination of strike-slip and thrust motions. Near the Caspian, 
continental crust is overriding the oceanic crust of the southern Caspian Sea. In 
Armenia there is both thrusting and strike-slip motion representing the eastward 
escape of the block, along with some thrusting. Further north on the greater 
Caucasus and northeastern shores of the Black Sea, earthquake mechanisms are 
almost purely thrust, as shown in Figure 9 (Alptekin et al., 1986). Taking into 
account the one thrust earthquake to the south along with the steep slope of the 
continental margin, it is possible that the eastern Black Sea is being overthrust by 
the continental crust, both from the north and the south. 

4. Mapping of active faults from remote Thematic Mapper and radar imagery (Gulen 
et al., 1991). 

While we have easy access to field sites in Turkey, access to northwest Iran and 
Armenia is more difficult. However, we have made considerable progress 
identifying active faults using space imagery (SPOT, Thematic Mapper and radar— 
SEASAT, SIR A, SIRB). 

GPS Densification of SLR Measurements in the Eastern Mediterranean 


The Turkey GPS experiment is pan of a broader effort, involving a number of 
institutions, to understand better the kinematics and dynamics of on-going tectonic 
activity in the Eastern Mediterranean region. The tectonics of this region are 
dominated by collision of the Arabian and African plates with the Eurasian plate. 
These plate interactions give rise to active subduction along the Hellenic arc, back- 
arc spreading in the Aegean sea and the Aegean trough region of S.W. Turkey, 
continent-continent collision and associated mountain building in the Caucasus- 
Bitlis Zagros ranges of eastern Turkey and adjacent areas, continental escape and 
major block rotation of the Anatolian and East Anatolian micro-plates with 
associated strike-slip faulting (including the North Anatolian fault), and strike-slip 
faulting along the Dead Sea transform. Virtually all types of continental deformation 
occur within the borders of Turkey. 

The MIT group, with support from UNAVCO and in collaboration with Hacettepe 
University, initiated GPS measurements in Turkey in 1988 with the establishment of 
14 sites primarily in western and central Turkey. During this campaign, GPS 
measurements were also made by our group at the 4 mobile SLR stations 



established in Turkey by WEGENER. In 1989, MIT/UNAVCO, in collaboration 
with the Turkish Union of Geodesy and Geophysics (TUJJB), established 18 GPS 
sites primarily in the continental collision-escape zone of eastern Turkey. 
Independent, but closely coordinated GPS campaigns were carried out 
simultaneously by IFAG/TUJJB, and the University of Durham/TUJJB. IFAG/TUJJB 
made measurements at SLR sites and established GPS “foot prints” on the scale of 
a few 10s of kilometers around each SLR site (~3 per SLR site). The 
Durham/TUJJB group observed a dense network in the Aegean trough region of 
western Turkey, including 5 of the sites established by our group in 1988. 

During the fall of 1990, regional GPS surveys to monitor tectonic deformation were 
undertaken by MIT/UNAVCO, IFAG/WEGENER, and ETH-Zurich. Excellent 
communication and cooperation were maintained throughout the planning and 
execution phases of the 1990 GPS surveys. Each group worked in collaboration 
with the TUJJB which also served to coordinate activities between groups. 
IFAG/TUJJB undertook a major campaign involving 10 receivers observing for 
approximately 6 weeks. The campaign included two phases: the first consisting of a 
regional campaign to observe 22 sites in the Turkish first-order geodetic control 
network and a number of points along a precise leveling line in western Turkey, the 
second consisting of a dense survey along a 250 km segment of the North Anatolian 
fault in western Turkey (east end of Marmara Sea to Gerede). ETH-Zurich, in 
collaboration with Istanbul Technical University, established a network of sites in 
the Marmara Sea region (west of IFAG North Anatolian fault network) to 
investigate activity along the westernmost segment of the NAF where it splays 
into a number of active strands. Details of the IFAG and ETH-Zurich experiments 
can be obtained from Hermann Seeger, IFAG and H. -G. Kahle, ETH-Zurich. 

Details of the M.I.T./TUJJB 1990 GPS campaign were reported in a previous annual 
technical report. Briefly, we reobserved 10 GPS sites in western Turkey 
established by our group in 1988. Simultaneous observations were carried out by 
IFAG at 3 SLR sites in central Turkey. In addition, reference marks on the scale of 
100 m were set and surveyed at 13 sites in western Turkey. Initial estimates of 
horizontal crustal motions determined from the 1988-1990 GPS observations are 
presented and described in Appendix 1. 

A first repeat of the eastern Turkey GPS network was accomplished in 1991. The 
principal objectives/accomplishments of the 1991 E. Turkey GPS Campaign included: 

1) Reobservation of the 1989 GPS network to obtain initial estimates of 

deformation rates and styles in the E. Turkey continental collision 
zone. 

2) Measurement of baselines crossing the western and central segments of 

the N. Anatolian fault to obtain additional information on rates of 
strain accumulation. 

3) Establishment of new stations to define better movements across the 

Dead Sea fault, and N-S shortening between the Arabian and 
Eurasian plates. 



4) Tying the E. Turkey network to the GPS network across the Caucasus 

established in July 1991. 

5) Establishing and surveying reference marks for E. Turkey GPS stations. 

6) Making initial GPS measurements at the new SLR pad in Karacadag (to 

replace Diyarbakir site) and tying the permanent Ankara GPS station 
to the Hacettepe University, Ankara site established by MIT in 1988. 

The 1991 field campaign was a collaborative effort with the TUJJB and was 
coordinated with GPS measurements carried out by IFAG. Seven Trimble 4000 
SST GPS receivers and 1 set of conventional survey equipment (for reference 
marks) were provided by UNAVCO for this experiment. Only 2 U.S. participants 
were involved in field operations (Reilinger, UNAVCO Engineer). All other 
participants (6 primary GPS operators, 8 assistants) were provided by TUJJB. 

Figure 10 shows regional GPS sites in Turkey established by our group between 
1988 and 1991. 

The 1988 and 1990 GPS observations in western Turkey (see previous section) 
were reduced at M.I.T. jointly with our Turkish collaborators. An abstract 
describing our initial results was submitted for presentation at the 1991 Fall AGU 
Meeting (Oral et al., 1991). The GAMIT and GLOBK GPS software packages 
were used to analyze double difference data to estimate site coordinates, orbital 
elements, atmospheric parameters, earth orientation parameters, and phase 
ambiguities. SLR site positions have not been imposed in our initial solutions. 
Figure 1 1 shows relative horizontal site velocities assuming no motion at Yigilca, 
an SLR station located north of the North Anatolian fault. The assumption of 
stability of the Yigilca SLR station is consistent with results reported from analysis 
of European SLR observations (Ambrosius, personal communication, 1991). The 
solutions shown in Figure 1 1 indicate right-lateral strike slip motion across the 
western and central sections of the North Anatolian fault at a rate of approximately 
2 1 0.5 cm/yr and southwestward “escape” of western Anatolia at approximately 4 
± 1 cm/yr. These initial results represent the first quantitative estimates of present- 
day deformation rates in this area. 
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Figure 1. The major deformational units of the Aegea, eastern Mediterranean. 

and Middle East. Large arrows show approximate movements of the 
Arabian and African plates with respect to the Eurasian plate. 







Figure d. Distribution of earthquake epicenters shallower than 50 km in the 
eastern Mediterranean and Middle East for the period 1961-1900 
(after Jackson and McKenzie, 1984). 




Figure 3. Major tectonic elements of Turkey (from Barka and Kadinsky-Cade, 
1988). 






Figure 4. Distribution of earthquake epicenters in Turkey for the period 1900 
1975. 
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Figure 6. Grid used in finite element models of Eastern Mediterranean region. 
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Figure 8 . Earthouake mechanisms in central and eastern Anatolia 

western Iran, and Caucasia. 




Figure 9 . Source Mechanisms of major earthquakes around the 

Black Sea. Note the thrust mechanism both on the 
northeast and southern margins. '(After Alptekin et al 
1986 ) 










Figure 11. Displacements (1988-1990) in western Turkey from GPS. 
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Abstract 


We present the first results of GPS measurements of regional crustal deformation in Turkey. 
The 1988 and 1990 GPS data from western Turkey were processed using the GAMIT and 
GLOBK packages. Derived "two-epoch- velocities" show a ~ 25 ± 9 mm/yr risht lateral 
motion across the North Anatolian Fault (NAF). Western Turkey moves southwest relative 
to the Eurasian plate at a rate of approximately 50 ± 20 mm/yr. Deformation models 
where piate boundaries are exoressea as contacts show rhat forces created by the collision 
of the African and Arabian piates with Eurasia are insufficient to explain the high rates of 
deformation indicated by GPS observations in western Turkey and SLR observations along 
the Hellenic arc. Southward migration of the Hellenic arc due to either "trench suction" or 
mantle convection induced by the down-going slab may account for the extra force driving 
the fast extrusion of southwestern Turkey. The relatively high elevation of eastern/ central 
Turkey may also contribute to the rapid extrusion of western Turkey. 

Introduction 

Deformation in the Eastern Mediterranean is controlled by the interaction of three plates — 
the African. Arabian and Eurasian plates. The bulk of the deformation in eastern Turkey 
is caused by the differential motion between the African and Arabian plates astainst the 
Eurasian plate, which manifests itself in the Bitlis and Caucasus convergent zones and crustal 
escape along strike-slip faults. This picture is compounded by extensional regimes in the 


9 



Aegean and western Turkey, subduction at the Hellenic and Cyprean arcs, and possible 
subduction in the Black Sea. 

SLR station coverage in the Eastern Mediterranean provides information on broad-scale 
plate motions, but is not sufficiently dense to delineate the details of deformation associated 
with ongoing continental collision. A GPS geodetic network has been established in Turkey 
to determine the rates and directions of relative plate motions and the internal deformation 
ot the Anatoiian plate. The collected phase and group delay data from the 1938 and 1990 
occupations in western Turkey were subjected to preliminary analyses using the GAMIT 
processing and the GLOBK adjustment packages: “two-epoch-velocities" were obtained. Our 
initial results provide the first quantitative estimates of present-day right-lateral strike slip 
motion across the central segment of the North Anatolian Fault and the southwest escape 
of western Turkev. 

Incorporating plate boundaries as contacts, we develop a finite element modelins scheme 
that facilitates the interpretation ot GPS and SLR measurements of crustal motions and 
provides a basis for investigating the dynamics of this complex region. These models demon- 
strate that N-S shortening across the collision zone is not sufficient to explain the observed 
deformations in western Turkey. Southward migration of the Hellenic trench and the ex- 
trusion of the Anatolian block are required to match Nuvel-1 plate velocities. SLR derived 
velocities at the Hellenic arc. and GPS site velocities. 
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Tectonic Background 


Severai authors have reported on the various geological and geophysical aspects of the East- 
ern Mediterranean where a number of coexisting tectonic regimes are confined ( Figure 1): a 
continent-continent collision zone (Bitlis-Zagros. Caucasus), an extensional province (western 
Turkey, northern Aegean Sea. southern Greece), major transform faults (the North. South- 
east and Northeast Anatolian Faults. Dead Sea Fault) and different staees of subduction 
(the Hellenic Arc. Black Seal. 

§engor and Yilmaz (1981). Robertson and Dixon (1984). $engor (1987), and Gealev 
(1988), provide a review of the geological evolution of the Eastern Mediterranean within 
the Alpine- Himalayan orogenic belt ana of its collisional development. The seismicitv, deep 
structure and active tectonics of the region are summarized in a number of studies (e.g., 
McKenzie. (1970. 1972); Jongsma 1975; McKenzie. 1976; Morelli. 1978: Papazachos and 
Comminakis 1978: McKenzie. 197S: Lort. 1978: §engbr. 1979: Camtez and Toksoz. 1980: 
Lvberis. 1984: Makris and Stobbe. 1984: Jackson and McKenzie. 1984; §engor et al . , 1985; 
Hempton. 1987; Barka and Kadinskv-Cade, 1988; Jackson and McKenzie. 1988; Philip et 
al.. 1989; Westawav, 1990 ). 

Present-day tectonic deformation (neo- tectonics) of the region is closely related to the 
northerly motion of the Afro- Arabian plate and medial to late Eocene events in the Red 
Sea (Hempton. 1987). During this period Africa and Arabia formed a single plate and the 
continental margin ol the Arabian platform collided with the Eurasian plate, closing the 
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back arc basin of the Tethvs. Thus, suturing began as subduction of Africa beneath Eurasia 
continued. As this process progressed into the early Miocene, the extension in the Red Sea 
and Gulf of Aden was accommodated by the separation of Arabia from Africa, resulting in the 
rapid convergence of the Arabian plate against Eurasia. Contemporaneously, the Dead Sea 
Fault initiated. By medial Miocene, the wide continental margin of Arabia was consumed, 
arresting its independent northerly motion. Subsequently, the first phase of extension in the 
Red Sea came to an end. as did the ieft-laterai motion along the Dead Sea Fault. Arabia 
and Africa, however, continued moving northward as a single plate causing excessive crustal 
thickening and shortening, and upiift in the convergent zone, until early Pliocene. The 
resulting convergent stresses were alleviated by the extrusion of continental wedges along 
transform faults (McKenzie, 19/2) . Therefore, between late Miocene and early Pliocene. 
Anatolian strike slip faults (the North. Northeast, and Southeast Anatolian Faults) were 
initiated (§engor, 19/9: Barka and Hancock. 1984) and the Anatolian piate started to wedge 
out of the zone of maximum compression. Such escape tectonics along transforms triggered 
Arabia to move independent of Africa. Consequently, the second stage of the Red Sea's 
opening commenced, providing a faster convergence rate for the Arabian plate. 

Global models for present-day plate motions using spreading rates, earthquake slip rate 
vectors, and transform fault azimuths show that the Arabian plate converges against Eurasia 
faster than the African plate (Minster and Jordan. 1978: Chase. 1978: Gordon and Jurdv. 
1986; De Mets et al., 1990). The differential motion between Arabia and Africa with respect 



to stable Eurasia, taken up along the Dead Sea Fault, amounts to 15-17 mm/vr I from 
De Mets et al.. 1990). The Alrican plate, which has been continuously movins northward 
relative to Eurasia since 95 Ma. converges on Eurasia at 11 mm/yr ( from De Mets et al.. 
1990). The Nuvel-1 data for the relative motion of Eurasia and Africa are obtained from the 
Azore- Gibraltar ridge, and that for Arabia and Africa come from the Sheba ridge. Because 
there are no data for the relative motion between Eurasia and Arabia, their relative motion 
is indirectly calculated from piate circuits. The data for Eurasia- Africa relative motion are 
from outside the Mediterranean, ^ince the Aubia-Somaiia motion is neglected, the spreading 
rates from the Red Sea do not constrain the Arabia-Africa motion. Further problems arise 
with the diffuse nature of regional continental deformation where the rigid plate assumption 
becomes less valid (McKenzie. 1977). The indirect methods of estimating plate motions in 
the Eastern Mediterranean and their probable uncertainties are important reasons for using 
SLR and GPS to obtain direct measurements of plate motions. 

The Arabian piate convergence against Eurasia created two tectonic flakes — the Van 
block (Iran plate) and the Anatolian plate (Figure 1). Bounded by the Caucasus and the 
Bitlis-Zagros sutures, and by the North. Northeast and Southeast Anatolian transforms, the 
Van block escapes eastward. The Anatolian plate continues to accommodate shortening 
by westward drift along the right-lateral North and left-lateral Southeast Anatolian Faults. 
The estimated rate of slip across the North Anatolian Fault from earthquake and geologic 
studies is quite uncertain, with values reported between 5 and 30 mm/vr (e.g.. ijengor. 1979: 
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Jackson ana Mckenzie. 1984: Barka and Kadinskv-Cade. 1988). West of the Mudurnu Valley 
the North Anatolian Fault gains an extensionai component I Barka and Kadinsky-Cade. 1988) 
and this feature becomes dominant in the North Aegean Trough (Lyberis. 1984). Although 
the central province of the Anatolian piate manifests relatively little internal deformation, 
the western province (western Turkey) is rich with extensionai features. The Buviik and 
Kii^iik Menderes Breakaway Zones, their cross horsts and grabens ( §engor et al.. 1985) 
demonstrate N-S and E-W extension, respectively. 

The bulk of the deformation in the region appears to have originated from the coilisionai 
process, i.e.. the forces applied at the plate boundaries. However, the driving forces for the 
present-day escape, as speculated by §engor et al. (1985) may also include forces caused by 
gravitational potential and/or basal traction. 

Space Geodetic Measurements 

SLR measurements 

The WEGENER-MEDLAS (Working Group of European Geoscientists for the Establish- 
ment of Networks for Earthquake Research Mediterranean Laser Ranging) project (Rein- 
hart et al.. 1985), which began in 1985, was the first step toward providing geodetic data 
that quantifies the deformations in this very active region. The first observational campaign 
started in 1986. Four sites in Turkey were visited during the third campaign: one site north 
of the North Anatolian Fault (Yigilca), two sites on the Anatolian plate (Yozgat and Me- 
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lengigiik). and another on the Arabian plate (Divarbakiri (see Figure ii. In 1988 and 1990 
these sites were occupied with mobile LRS (Laser Ranging System! and laser ranee measure- 
ments to LAGEOS were performed. Analyses of these normal point data ( Sellers et ai. 1991: 
Noomen et ai. 1991: Wakker et ai. 1991: Cenci et ai, 1991) show a 20 ±10 mm/yr right 
lateral motion across the North Anatolian Fault and a 50 ± 20 mm/yr southwest motion of 
Aegea relative to Eurasia. 

GPS measurements 

In 1988 with support trom UNAVCO and in collaboration with Hacettepe University, 

initiated GPS measurements in Turkey to establish fourteen sites primarily in western and 
central Turkey (Figure 1). During this campaign, GPS measurements were also made at the 
four mobile SLR stations established in Turkey by WEGENER-MEDLAS. 

During the summer/ fall of 1990. regional GPS surveys were undertaken by M.I.T.. IFAG. 
and ETH-Zurich in collaboration with Turkish Union of Geodesy and Geophvsics (TUJJB). 
During this campaign the M.I.T. group reobserved ten GPS sites in western Turkey that had 
been established in 1988 (Figure 1). Simultaneous observations were carried out by IFAG at 
three SLR sites in central Turkey. 

GPS campaigns were also undertaken in 1989 and 1991 by M.I.T. and Turkish collabo- 
rators which focused on deformation in the collision zone of eastern Turkey and neighboring 
areas in the Trans-Caucasus (Reilinger et al.. 1991) . Figure 2 shows all GPS sites established 
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by the M.I.T. group tor the period 1088-1991. Data from the eastern Turkey campaigns are 
not included in this study. 

Data 

Figure 1 shows GPS sites common to the 1988 and 1990 campaigns. Baselines in the west- 
ern Turkey network range from 95 km to 1200 km. with an average baseline length of 350 
km. In both campaigns, observations were made at thirteen sites: Akdag, Ankara. Antalya. 
Cesmeihca. Cine, Kandilli. Ismetpa§a. Mekece. Melengigiik. Pamukkale. I'ludag, Yigilcaand 
Yozgat. No repeat observations at the Avvahk. Demirkoy. Divarbaktr. Ka§ and Odemi§ sites 
were taken in 1990, although their occupation is planned for 1992. All receivers and antennas 
used for the 1988 survey were TI-4100s. The repeat experiment used a mixture of TI4100, 
MiniMac. and Trimble 4000s (SST and SLD), introducing asynchronous and irregular sam- 
pling of phase and group delay. During the experiments, selective availability was off and 
anti-spoofing was not introduced. In 1988 oniv six satellites were available (prns 3. 6. 9. 11. 
12 and 13). The repeat experiment observed four additional new satellites iprns 2. 16. 18 
and 20). 

GPS Data Processing 

The GAMIT software package (King and Bock. 1991; King et al., 1985) was used to process 
the 1988 and 1990 GPS data. GAMIT is a suite of programs that fulfills a number of 
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tasks, such as: (1) converts the field data into internal formats: (2) generates tabulated 
reference orbits and their partials for satellites irom broadcast / adjusted ephemeriaes with a 
celestial model: (3) calculates the theoretical observations based on a geometrical model: (4) 
determines partial derivatives and differences from observed phase measurements: (5) detects 
outliers and gaps, ana edits one-way data: and (6) estimates three-dimensional relative 
positions of observational sites and orbits, and atmospheric clock and ambiguitv parameters 
irom a weighted least squares analysis ot phase delay measurements at the GPS frequencies. 
A simple flow cnart of GAMIT soitware is uiven in Figure 3. 

Solutions are carried out in the SVo reference frame (Murray, 1991) derived from VLBI 
(GLB659) and SLR (CSR8902) with the AM 1-2 hot spot model (RM2). Carrier beat phase 
observables and pseudoranges are used, and geometric terms are modeled (King et al.. 1985: 
Feigl, 1991). Sea level conditions are applied to determine the atmospheric delav at the 
zenith. The Saastamoinen 1 19(2) model is applied for dry and wet zenith delav. The mapping 
function for delays is definea by Davis et al. (1985). Ionospheric delay is not included in the 
model, however it is calculated explicitly at the solution stage by estimation or application 
of the LC observable, or a combination of both as an ionospheric constraint. 

Satellite clock parameters (offset, rate and acceleration) are modeled with a quadratic 
polynomial. No relativistic effects are included. Receiver clock parameters (offset, rate, 
acceleration and rate of acceleration ) are estimated from pseudoranges epoch by epoch using 
a cubic polynomial with jumps. 
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Merit standards with an updated GM value and IERS tables are used for orbits and 
expressed in an inertial coordinate system. Perturbed satellite motion is defined with six 
osculating elements (cartesian position and velocity). The earth’s gravity field is determined 
by the GEM-L2 model up to degree and order 3. The gravitational attraction of the sun and 
moon are considered. Non-gravitationai parameters acting on the satellites are absorbed in 
the direct radiation pressure, x-axis and v-axis biases. A tabular ephemeris is integrated for 
each session at intervals oi 22 minutes. For a single session the short-arc solution orbital 
trajectory is ~ 12 hours long. For a long-arc single session '"global'’ solution. 3-4 day long 
orbital trajectories are integrated and applied to a single day solution. 

Removal of outliers and cycle-slip repairs are performed automatically with one-way, 
between-satellite difference, and using double difference data. Manual editing under X- 
windows is routinely practiced to eliminate cycle-slips overlooked during the automatic clean- 
ing procedure. 

A weighted least-square scheme is used to estimate the delays in the propagation medium, 
clock and ambiguity parameters (wide-lane [L2-L1] and narrow-lane [Ll] biases;, 3-D site 
positions, and satellite state vectors (six Keplerian elements, three non-gravitational param- 
eters). Each single day solution consists ol four steps (Feigl, 1991). First, all ambiguity 
parameters are set free and the ionosphere-free observable LC is used to estimate all the 
parameters i Bock et al., 1986). Second, keeping all the parameters other than those whose 
ambiguities are fixed to the values obtained in the earlier step, wide-lane biases are estimated 


11 



to fix the ambiguity parameters to integers. LI and L2 phases, and wide-lane combinations 
(if and when pseudoranges are available) are used in this step. Third, narrow lane biases 
are estimated. The last step involves solving for all the parameters except wide-lane and 
narrow-lane biases which were successfully fixed to an integer. This is a ••biases-fixed' ! solu- 
tion. Biases-fixed and biases-free solutions are repeated for sufficiently loose constraints to 
obtain estimates and variance-covariance matrices as an input to the "global" solution. 

A global solution requires estimates and variance-covariance matrices of lons-arc single 
session solutions. The GLOBK software package (Herring et al.. 1990: Herrins, 1991), a 
Kalman filter program, combines the variance-covariance matrices of these solutions to es- 
timate the 3-D site positions and velocities, satellite state vectors, and earth orientation 
parameters (polar motion and UTl-UTC). 

Application to Data 

The Onsoia. Richmond. Wettzeii. Westiord. Tromso. and Yellowknife CIGNET sites were 
used as fiducials in processing the 1988 campaign. The repeat campaign lacks Yellowknife 
anc ^ Tromso data due to data outages. Single day short-arc solutions have approximately 12 
hour long orbits. Though the fiducial concept has been followed, instead of fixing fiducial 
site positions, we applied tight constraints on their positions (1-10 mm). Global solutions 
are performed with 3-4 day long arcs. Lsing the GLOBK package, precise determination 
of satellite state vectors, adjustment of site coordinates, and estimation of velocities are 
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attempted. 


The solutions presented are not **biases-iixed : ‘ and may still suffer from inhomogeneous 
fiducial site distribution and unconstrained orbital parameters. Repeatability in the east, 
north, and vertical is 10-30 mm + 1 parts in 10 8 . 


Preliminary Results of GPS Data Processing 

To understand the interaction between the Eurasian. African and Arabian plates, site veloc- 
ities expressed relative to the Eurasian plate provide a favorable frame of reference that can 
be conveniently described by the Onsola, Tromso and Wettzell sites. The baseline lengths 
between these and local sites are on the order of a few thousand kilometers. Thus, their 
relative errors, which translate to decimeters, do not allow the local sites ? velocities to be 
resolved. Therefore, we adopt the frame defined by the Yigilca site which expresses “two- 
epoch-velocities relative to the northern block of the NAF. This approach can be tested by 
examining the relative motions at Ismetpa§a. a site 25 km east of Yigilca. also located north 
of the NAF. Since there is no significant motion at Ismetpa§a at 95% confidence intervals, 
we conclude that rates obtained relative to Yigilca represent an adequate frame of reference. 
Furthermore, SLR solutions indicate no significant motion of Yigilca relative to stations in 
Europe (Ambrosius et a/., 1991), indicating that within the resolution of our measurements. 
Y lgilca is a part of the Eurasian plate. Thus, the velocities derived in this study are in a 
fixed Eurasian plate reference frame. 
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Table i provides geocentric site coordinates l WGSS4 datum), horizontal ana vertical rates 
and their standard deviations and correlations. Figure 4 shows the "two-epoch-velocities" 
relative to Yigilca. The 95% conridence levels are high and vary between 10 and 20 mm/yr. 
Subsequently, rates are unresolvable at some sites due to the bias-free nature of the solution, 
residual orbital errors, and inadequate and non- homogeneous fiducial coverage. 

The sense of motion indicated by the GPS rate vectors is generally consistent with those 
determined from geological observations. The magnitude and direction of GPS rate vectors 
are also consistent with some of those aerived from SLR observations in the Aegean within 
their 95% confidence intervals, although some discrepancies remain between SLR solutions 
reported by various groups (e.g.. Sellers et ai. 1991: Wakker et al., 1991). 

Right-lateral motion across the North Anatolian Fault is easily detectable, with the 
preliminary GPS estimate across this transform being 25 ± 9 mm/yr. The Anatolian block 
is rapidly escaping SW at a rate of 50 ± 20 mm/vr. The rates across the Marmara Sea 
demonstrate that the right-lateral motion continues westward. The significant motion at 
Kandilli. relative to Yigilca. may reflect the extensional character that NAF assumes west of 
the Mudurnu Valley. Higher velocities at Mekece and Uludag also indicate that the motion 
across the southern branch of the western segment of the NAF is different from that on the 
central segment. Though relatively high confidence levels make it hard to resolve the N-S 
extension in western Turkey, there is at least a 10 mm/yr motion. A better estimate may be 
reached after the planned 1992 repeat measurement at the Odemi§ site on the Bozdag Horst. 
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Modeling and Discussion 

A finite-eiement modeling scheme, constrained by the SLR and GPS measurements of crustal 
motions as well as geological and seismological observations, was developed to investigate 
the kinematics and dynamics of the interaction between the Arabian. African, and Eurasian 
plates. The equation of equilibrium is solved for an aggregate of 2-D elastic plates (plane 
stress I in contact. A general potential energy equation in variational form is derived for fric- 
tionai/ tensional contacts lOral t 1 aL. i991L It includes body forces, basal traction, ambient 
strain/ stress, and variable piate thickness and is expandable to 3-D. This non-linear contact 
problem is iteratively solved for present-day plate velocities and forces arising at the contacts 
(Wang and Voigt, 1969). Similar contact algorithms (Bathe and Chaudhrv, 1985: Chaudhry 
and Bathe. 1986; Bathe and Mijailovich. 1987) and joint element approaches (Goodman et 
aL ' 1968: Goodman. 1975) are. in a broad sense, equivalent. Contacts are modeled with dual 
nodes, one node on each plate ( Kasapoglu and Toksoz. i9S3). Compatibility is assured by 
Coulomb-Mohr failure criteria defined by three modes of contact allowing sudden velocity 
jumps across predefined boundaries sticking (continuum), sliding (slip), and tension release 
(separation). The velocity field is discretized by 4-node quadrilateral serendipitv elements. 
Gaussian quadrature integration is carried out at 3 x 3 Gaussian points. 

Forces other than those at the contacts are ignored in the current models. Major plate 
boundaries are modeled as contacts: Dead Sea. North. Northeast, and Southeast Anatolian 
Faults. Aegean Trough, Aegean Block Faulting, Mediterranean Subduction. Arabian Colli- 
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sion Zone and Greater Caucasus Thrusting ( Figure 5ai . Due to the 2-D nature of the model, 
subduction can only be modeled explicitly such that the convergence across the Hellenic arc 
has to be specified. The contacts, assumed to be very weak, have almost no cohesion and 
have a low friction coefficient. 

We have considered several deformationai models on which Nuvel-1 velocities for Arabia 
and Africa relative to fixed Eurasia are imposed. The African and Arabian plates are assumed 
to be rigid with uniform velocities relative to the Eurasian plate of 11 mm/vr fat X 30. 
E 45) and 27 mm/vr (at X 30. E 35). respectively. These forcing conditions, applied on and 
within these plate boundaries define, by default, the deformations caused by the collision 
of Arabia against Eurasia, provided that the Anatolian block is treated as a separate plate 
and the Hellenic arc is explicitly modeled. The velocity field obtained from such a model 
underestimates the rate of deformation across the NAF and in western Turkev. Instead, we 
combine velocities obtained from GPS and SLR observations and model the Hellenic arc 
moving WSW at a rate of 60 mm/ vr (30 mm/yr east of Kattavia SLR site. Rhodes), relative 
to Eurasia. This implies the application of some boundary forces on the overriding plate. 
Figures 5b and 5c show the velocity and strain fields, respectively for a model with Nuvel-1 
plate motion and prescribed SW motion of the Hellenic arc. Sidewavs escape is remarkably 
dominant. Two contrasting strain regimes are noticed — extension in western Turkey and 
shortening in eastern Turkey. With an average strain rate of 8 x 10 -6 yr _1 . the Van block 
escapes XE at ~ 15 mm/yr. The motion across the NAF and SEAF is predicted to be ~ 20 
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mm/vr. The contact along the NAF does not yield a pure slip mode. There is a component 
of tension release. The deformation across the NEAF is less than iO mm/yr. Our models 
also predict an opening at the Karliova and Maras triple junctions, and in the Marmara 
Sea. The extensional regimes in the Aegean Sea and southern Greece are consistent with 
geological observations. 

The velocity field is slightly lower than that of GPS. suggesting the necessitv of modeling 
a region inclusive of the Aegean Sea. as well as the Hellenic arc. to produce GPS velocities 
in western Turkey. This implies that basal slip should be introduced to the models. Another 
option is to increase velocities at the Hellenic arc. This is done, however, at the expense of 
exceeding SLR determined velocities. 

Comparison of Figures 4 and 5b indicates that the velocity at Melengiflik deviates sig- 
nificantly from that predicted by our model. This may be due to the motion of the block 
bounded by the SEAF and the Teke wedge (the wedge is created as the Mediterranean 
subduction jumps from the Hellenic arc to the Cyprean arc). The boundaries of this block 
and wegde are not included in our models. Instead, a continuous boundary is assumed 
(Figure 5b). Although this is a deficiency in the current model, it shows how the contact 
approach, by allowing velocity jumps across boundaries, can correctly define the plate mo- 
tions. A similar argument is valid for the breakaway zones and their cross horsts and grabens. 
whose features must be modeled as contacts in order to delineate the relative motion across 
each small block. 
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Conclusions 


The GPS measurements in Turkey, combined with more widely spaced SLR coverage, provide 
valuable information for quantifying the diffuse deformation in the Eastern Mediterranean. 
Site velocities obtained by GPS agree with those of SLR within their 95% confidence levels. 
The accuracy of the velocities obtained from the 1988 and 1990 campaigns will be improved 
as the 1989 and planned 1992 campaign data are included in the anaivsis. 

Site velocities obtained irom two epoch measurements show a rapid escape of western 
Turkey along the Anatolian transiorm rauits (50 ± 20 mm/yrj. The relative motion found 
across the central segment of the XAF (25 ± 9 mm/yr) provides the first quantitative 
estimate of present-day strain accumulation with implications for seismic risk assessment in 
this heavily populated region. Likewise, the geodetic data show that the western segment 
of the XAF, splintered in several branches, is much more active than was expected. The 
rates obtained at Mekece mav retlect the motions of the Gernlik block, which is bounded by 
these branches. The extensive deformation in and around Bozdag Horst, quantified by GPS 
measurements, demonstrates the active extensional tectonics of this region. 

Numerical models show that the deformation in western Turkey cannot be accounted 
for by collisional processes alone. A corner flow model driving the Aegea. combined with 
gravity driven processes (YVdowinski et a/.. 1989), may explain the high rates of relative 
motion observed in western Turkey. It is shown that modeling plate boundaries as contacts 
across which sudden velocity jumps are allowed is an adequate technique to account for slip 


18 



and tensionai modes. The limitations posed by explicit modeling oi dip siip can be removed 
by expanding to the full 3-D. 
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Figure Captions 


Figure 1 : Tectonics of the Eastern Mediterranean region (after Sengor et al. . 1985. 
modifications by A. A. Barkai. Sites occupied during the 1988 and 1990 campaigns are 
indicated by solid circles. Open circles show sites occupied in 1988 only. 1: Ankara. 
2: Yozgat. 3: Melengiplik. 4: \ lgilca. 5: Ismetpa§a. 6: Mekece. 7: Kandilli. S: Demirkov. 
9: lludag, 10: Akdag. 11: Ayvahk. 12: Odemis. 13: (^ine. 14: Pamukkale. 15 : Ka§. 
16: Antalya. 1 . ; Diyamakir. .i: (,'esmemca. 

Figure 2 : MIT GPS network in Turkey and Trans-Caucasus. Sites occupied between 
1988 and 1991. 

Figure 3 : Flow chart of GAMIT software. 

Figure 4 : The two-epoch-velocities relative to Yigilca. Site velocities are obtained 
by combining variance-covariance matrices of long-arc single-day GAMIT solutions in 
GLOBK. Yigilca is denoted by a large dot. 

Figure 5a : The finite-element discretization and plate boundaries as contacts. Solid 
lines show contacts: Dead Sea. North. Northeast, and Southeast Anatolian Faults. 
Aegean Trough, Aegean Block Faulting, Mediterranean Subduction. Arabian Collision 
Zone and Greater Caucasus Thrusting. Scale unit is meters. 
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Figure 5b : The velocity field. Boundary conditions: Eurasian piate is held fixed: 
Arabia and Africa move 27 mm/vr and 11 mm/vr relative to Eurasia, respectively: 
Hellenic arc moving WSW at a rate of 60 mm/ vr (30 mm/ vr east of Kattavia. Rhodes), 
relative to Eurasia. Arrows are at the end of the velocity vectors. At the contacts, 
there are two velocity vectors for each plate. Scale unit is meters. 

Figure 5c : The strain field calculated at the Gaussian points of each [finite! element, 
inward directed arrows indicate shortening. Scale unit is vr -1 . 



Table 1 Horizontal and Vertical Rates i '‘two-epocn- velocities ’ ) relative to Yigilca. Site coordinates axe geocentric. 


Rates Rate 


SITE 

Longitude 

Latitude 

East 

North 


& nn 


Vertical 

-■* 


(deg) 

1 deg) 

(mm/yr) 

(mm/yr) 




(mm/yr) 


AKDACS 

27.873 

39.006 

-■55.0 

-13.2 

7.1 

4.0 

0.707 

-122.2 

10.2 

ANKARA 

32.726 

:o.87i 

-22.6 

-10.6 

? ° 

1.6 

0.581 

-132.8 

5.4 

ANTALYA 

30.583 

6.822 

-13.0 

-14.7 

7.1 

n -» 

0.585 

-137.0 

10.5 

QE§MEILICA 

26.385 

8.311 

-58.1 

-29.9 

7.3 

5.0 

0.632 

-142.9 

9.5 

CINE 

23.081 

* 609 

■41.1 

-27.0 

6.9 

4.3 

0.651 

-139.7 

10.5 

ISMETPA^A 

32.570 

• > 881 

'-.6 

0.7 

El 

1 . 7 

0.631 

- 1 23.8 

" . 2 

MEKECE 

30.026 

;0.465 

-48.7 

-10.0 

5 . 6 

2.9 

0.736 

-125.6 

10.6 

MELENGIQLIK 

33.191 

“.378 

S.l 

-13.7 

6.8 

2.3 

0.429 , 

-21.3 

10.9 

KANDILLI 

29.064 

11.064 

-13.7 

5.2 

4.4 

2.6 

0.705 

-99.7 

8.3 

PAMUKKALE 

29.136 

37.941 

-23.5 

-12.5 

7.0 

3.8 

0.675 

-144.5 

10.5 

uludaG 

29.141 

50.122 

-35.5 

0.2 

5.8 

2.9 

0.754 

-150.3 

8.8 

YOZGAT 

34.813 

39.801 

-10.2 

-4.8 

6.9 

3.3 

0.583 

-108.4 

10.0 
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